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Abstract 


The  existence  of  an  additional  symmetry  element 


allows  one  to  reduce  the  size  matrix  needed  to  cal¬ 


culate  the  vibrational  normal  modes.  In  this  paper 
we  detail  how  to  reduce  the  matrix  appropriate  to  an 
alternating  DNA  polymer  which  has  a  symmetry  about 
the  x-axis. 


I.  Introduction 

y*  A  considerable  number  of  observations  both  by  infrared  absorption 

^  t -  ■  jr  t  ^  ^  <  '  A  ^ 

and  Raman  scattering  have  been  made  on  samples  of  double  helical^t)NA^ 

Many  of  these  observations  are  on  DNA  polymers  which  have  or  come  close 
to  having  structural  symmetries.  We  have  carried  out  helix  lattice  normal 
mode  calculations  for  a  number  of  repeating  DNA  polymers.^  The  agreement 
between  the  theoretical  predictions  and  experimental  observations  has 
been  quite  good  and  a  reasonably  good  set  of  force  has  been  refined 
for  such  calculations.  Most  such  calculations  have  been  on  the  simplest 
repeating  homopolymers There  are  a  considerable  number  of  observations  on 

(v- 

alternating  homopolymers  and  we  have  also  calculated  the  helix  lattice  modes 

2 

for  a  number  of  these  cases.  The  alternating  polymers  have  a  unit  ceil 
consisting  of  two  base  pairs  which  is  twice  as  big  as  the  simplest  homo¬ 
polymers.  A  true  helix  can  be  generated  by  applying  a  screw  axis  operation 
to  the  unit  cell  many  times. 

In  applying  the  helix  lattice  formalism  one  gets  the  normal  modes  by 
diagonalizing  a  matrix  which  has  dimension  of  3  times  the  number  of  atoms 
in  the  unit  cell.  The  matrix  for  alternating  DNA  is  thus  a  square  array 
twice  as  big  in  both  columns  and  rows.  The  structure  of  the  double  helix 

is  such  that  the  alternating  DNA  has  an  additional  symmetry,  a  C 2  operation  ,  

or 

that  Is  not  an  exact  symmetry  in  the  simpler  homopolymers.  The  point  of 
this  paper  is  to  detail  how  one  can  make  use  of  this  additional  symmetry  to 

>n. 

reduce  the  size  matrix  needed  to  be  diagonalized  for  the  alternating  homo- 


po 1 ymers . 


_D 1 a t  r i bu  1 1 on/ _ 

Availability  C«d«f 
!Av«iL«n#/*r~ 
tu-t  s»«ni  \ 

_  +  \  (  ) 

At  \mhi  f‘CTiV 


■ras 


i-.  t  .  r.  t  «.i  i.c  ».»a. ...  c. 


The  method  has  been  used  to  reduce  the  caiculat  ion  for  z-form 
po  I  y  (dC-dC)  •  |>o  I  y  (dt'.-dC)  and  il  has  led  in  a  reduction  in  machine  I  ime 
hy  about  an  order  of  magnitude  in  the  calculation  of  the  entire  normal 
mode  spectrum.  In  an  attempt  to  describe  the  method  for  incorporating 
the  additional  symmetry  in  as  general  a  manner  as  possible ,  we  do  not  use 
specific  details  of  the  internal  lorce  field  of  the  alternating  polymers. 
Many  details  specific  to  DNA  will  be  left  to  a  later  publication. 

The  paper  is  arranged  as  follows:  in  section  1 1  we  discuss  the 
general  features  of  helical  structure  and  establish  the  dynamical  matrix; 
Section  III  is  devoted  to  a  C.,  symmetry  discussion  and  section  IV  develops 
scattering  selection  rules;  conclusions  are  in  section  V. 
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II.  General  Properties  of  Heiicai  Structure 


Let  us  consider  an  infinite  coin  molecule  with  helical  structure 
such  that  each  unit  cell  is  transformed  geometrically  into  the  next  by 
the  operation  (R|p),  where  R  is  the  rotation  through  an  angle  i|/  about  the 
helical  axis  2: 


R  = 


cosp 

simp 


0 


-sin p  0 
cos'J/  0 
0  1 


(2.1) 


f  * 

and  p  =  pz  is  the  translation  along  z.  If  the  unit  ceil  is  the  smaiJest 
repeating  unit  of  the  molecule,!'  is  the  pitch  angle  and  p  is  the  pitch. 

The  structure  of  the  molecule  is  invariant  under  operation  (R  |np),  for 
if  X™  is  the  equilibrium  position  of  i^  atom  in  m^  unit  cell,  then 


(Rn|np)X^  =  RnX™  + 


npz 


►m+u 

Xi 


(2.2) 


The  choice  of  a  unit  cell  is  not  unique,  usually  one  takes  the  unit 

cell  as  the  smallest  repeating  unit  or  as  the  primitive  lattice  cell  of  the 

3 

corresponding  space  group.  With  the  choice  of  a  unit  ceil  ot  S  atoms  one 
gets  3S  vibrational  bands.  Each  band  has  vibrational  modes  that  differ  by 
having  a  phase  shift  from  one  unit  ceil  to  the  next.  The  number  of  phase 
shifts  is  equal  to  the  number  of  unit  cells  in  the  macroscopic  sample, 
approaching  a  continuum  of  phase  shifts  as  the  number  of  cells  in  the  sample 
approaches  infinity.  The  different  phase  shifts  are  equivalent  to  different 
waves  of  varying  wavelengths.  Any  modes  with  varying  phase  are  infrared 
and  Raman  inactive.  ^  In  cases  where  the  original  symmetry  of  the  system 


is  broken,  inactive  modes  mav  become  active. 


-4- 


In  our  helix  lattice  calculations  the  unit  cell  is  chosen  as  the 

1  ,2,  “5,6 

smallest  repeating  unit  as  tins  gives  rise  to  the  smallest  matrix 

for  solution.  In  the  harmonic  approximation  the  potential  energy  of 
the  molecule  is 


..  .  Mrr--mn  *n 

'  >■  q  .  •  D .  .  •  n  . 

.  .  i  1 1  i 
mil  l  ) 


(2.  J) 


where  m,n  label  the  unit  cells  and  i,j  label  the  atoms  in  a  unit  cell; 

q*1"  =  /m,  7  is  muss  weighted  cartesian  displacement  where  m  =  mass  of 

atom  i  and  u™  is  the  displacement  of  atom  i  in  m*1^  cell  from  its  equilibrium 

position  X™  is  the  dynamical  matrix  containing  the  force  constants 

nn 

between  the  atom  at  X^  and  the  atom  at  X..  The  kinetic  energy  is 


T  -  1  i  ^  ^ 

T  -  2  \  Vqi 

mi 


(2  4) 


And  the  equs.  of  motion  could  be  written  as 


,  cHim  mi  2  >m 
V  D .  .  -  q  .  =  ie  q  . 

M  l  i 

it  I 


(2.5) 


—  ioit 

whole  tin.'  harmonic  behavior  e  is  assumed. 

Bv  using  Born-von  K.irmnn  boundary  condition  and  taking  into  account 
the  helical  symmetry  we  impose  a  solution  which  is  a  mode  of  no  phase 
s  h  i  I  to: 


>m 

q  ;  = 


1 


m •  .  .  im  ' 
d"  R  q  ,  (u  )e 


(2.6) 


whe  1 1 


'!,(">  - 

1  A> 

y  Zn  m 


,  ,  -m  -  im1* 

/.  K  q  .  e 


(2.7) 


s.it  ist  ios  the  secular  oqn. 


whe  re 


)■  D.  .  *q  .  (•>)  =  .  (0) 

j  1.1  I  l 


5*.<e)  -  j:  R-"*o»Ve'u"“n)0 

lj  mu 


(2.8) 


(2.9) 


provided  that  R  ^  D^R11  is  only  the  function  of  m-n  for  a  certain  pair  of 
(i.j).  In  other  words,  R  10  ^o'j'^R11  is  unchanged  under  operation  ( E  |  p )  where  E 
is  identity.  The  reality  of  q™  implies 

q^-'O  =  q*(0)  ,  (2.10) 

i  — ► 

uni  i  hr*  reality  of  the  potential  energy  implies  D.  .  is  hermit  ian: 

4— ►  « -►  + 

Dlj=Dji  ■ 

For  the  molecule  possessing  no  symmetry  operation  which  takes  0  to  -0, 

eqn(2.10)  is  simply  the  statement  of  double  degeneracy  due  to  time  reversal. 

12  5  6 

It  has  been  pointed  out  ’  ’  ’  that  the  f ree  center  of  mass  motion  ol  the 

helical  molecule  will  occur  at  "  =  0,  +  iji.  The  eigenvectors  of  the  free 

motion  can  be  obtained  from  eqn.(2.7)  and  are  listed  in  table  !.  Free 

■enter  of  mass  motion  is  independent  of  force  constants  and  gives  restrictions 

on  d  v ii. inti  im  I  ilia  1 1  i  >,  0  : 

i  I 


>:  /m.  D  1  =  0,  0  =  0 

j  ttz 


C>.l  2.0 


>:  /m.  (D 1  j XJ  -  DlJx')  =  0,  0  =  0 

]  ay  x  ax  y 


(2.12b) 


>:  /m.  (D  J  +  i  D  1  )  =  0,  0  =  +  iji 

jjax—  ay  — 


(2.12c) 


1  I  1  - - 

whore  D^,  =  (D  , =  (X^)^  ,  and  a,B  =  x,y,z;  the  factor  /m^  comes 

I rom  the  fact  that  q^(fJ)  is  mass  weighted. 

8  -► 
From  n  classical  point  of  view  the  induced  electrical  dipole  moment  u 

of  the  molecule  is  given  by 


Mn  mu,  r  — 

|i  =  X  e.u.  =  >.  e.q./rm.  , 

li  .11  l 

m, l  m,i 


(2. J  J) 


where  o.  is  the  atomic  partial  charge  assigned  to  atom  i.  For  the 

-un 

tree  center  of  mass  motion,  tor  example  the  longitudinal  translation  u^  =  z, 

then  is  no  induced  dipole  and  eqn.(2.13)  implies  E  e.  =  0,  i.e.,  the 

1  1  10 

unit  cell  ot  the  molecule  must  be  neutral.  It  is  well-known  that  the 
dissolved  DNA  molecule  is  not  neutral.  To  allow  free  motion  of  DNA 
in  solution  the  counter  ions  must  be  included  in  the  DNA  to  make  the  unit 
,(.•11  neutral.  Otherwise  the  free  motion  of  DNA  in  solution  is  not  allowed. 


'.'.V.VA' 
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III.  The  Symmetry 


The  standard  z-form  DNA  poly(dG  -  dC)*(dG  -  dC)  has  a  symmetry  axis 
which  is  perpendicular  to  the  helix  axis  and  is  usually  taken  to  be  the  x 
axis  in  cartesion  coordinates.^  In  such  a  coordinate  system  an  atom  in  the 
m^'  coll  on  one  strand  has  equilibrium  position  (x,  y,  z),  then  the  corre¬ 
sponding  atom  in  -m^1  cell  on  the  antiparallel  strand  locates  at  (x,  -y,  -z) . 

For  those  DNA  molecules  without  C ^  symmetry  the  backbone  atoms  still  have 

„  1 2 
approximate  G  symmetry  and  x  is  picked  in  the  same  way 


For  this  symmetry  discussion  we  rewrite  the  secular  eqn.  (2.8)  as 


loll  ows 


.V  e-V  ->  2-* 

V  nii'qj  +  T;  Did,qS.  =  w  qi  ’ 

.1  ^ 

*W  1  Vv*J%  • 

It 


(3.  I) 


such  that  the  dynamical  matrix  D  is 


4  V  «  v 

D,  •  *>,  „ 

k  i  M 


(3.2) 


win.'  re  i,  i  -  I.....S  label  the  atoms  on  one  strand;  k  '=  i  +  S  ,  V  -  j  +  S  , 

a  a  .i 

aid  k,  t  label  the  items  on  the  antiparallel  st i  md  in  such  a  wav  that  atom  k 
is  the  same  atom  as  atom  i  but  on  a  different  str.ind.  S  is  the  number  of 

.i 


atoms  to  on 


e  strand.  3S  ~  S  number  of  atoms  in  i  unit  cell.  It  has  been 


pointed  out  hv  Lax*  that  the  operations  that  do  not  leave  invariant  provide 
relationships  between  the  corresponding  matrix  elements.  Tn  our  case  t  hi  C,, 


■i  at  ion  inverts  t  he  orientation  of  DNA  and  hence  takes  h  to  then  the 


to)  lowing,  relations  applv: 


(3.3) 


V'0  ■  *2*  R2x  *  '6kt<">"  R2x 


Vi.<“>  "  R2x  'Dk)<-'l)R2x  “  R2x\i(">‘  R2x  ' 


when* 


R2x  ■  R2x 


,ind  fK  .  (-H)  =  D^.(O)  follows  immediately  from  (2.9). 


F.qn.  (3.3)  implies  that  the  solution  of  secular  eqn.  (3.1)  takes  the 
I orm  of 


V°>  =  R2x  ^ 


(3.4) 


where  i>  is  a  common  arbitrary  phase  angle,  and  could  he  chosen  to  be  zero, 
Mich  that  the  real  part  of  eigenvector  is  symmetric  and  the  imaginary  part  is 
antisymmetric  under  C  operation.  By  imposing  the  solution  (3.4)  with  <p  -  0 


and  separating  Re  q  and  1m  q.  as  independent  variables  the  secular  eqn.  (3.1) 


becomes 


4  ->  -4  2~k  . 

J:  Vqi +  i:  (Vq»  - w  q*  * 


« 


(3.3) 


/v<i’ +  rsas  ■  ,"2<:k 


where 


1  I 


Re  |V.  .  +  R"1  V  R0  1  , 
ii  2x  l ?  2x 


-  I  m  f  I) .  -  R  \).  .  R  | 

l  •  2x  ii  2x 


kj 


-  «  V 


R(|"k.  -R2x  '>MR2XL 


'I  2 


*  ^  I  ^ 

Re  q,  ,  q  =  Re  q.. 


And  the  hormitian  matrix  D  (eqn.  ( 3 •  2))  becomes  a  real,  symmetric  matrix  C. 


KcfV.  .  +  R  1  V  R.,  ]  ,  -ImfD..  -  R"1  V.  .  R„  ] 

lj  2x  if  2x  if  2x  l  ]  2x 


‘■•Kj  +  R2x\,B2J  ’  Re(V  -  R2«%  V 


At  0  =  0,  I)  is  reaL  and  tlius  C  is  block  diagonal.  Eqn.  (3.5)  becomes  two 
independent  eqns,,  on<'  for  the  symmetric  and  one  for  the  antisymmetric  eigen¬ 
vectors.  By  checking  table  I  we  find  that  at  0=0  the  acoustic  longitudinal 

mode  and  torsional  acoustic  mode  are  both  antisymmetric  under  operation.  Tin- 
*  * 

restrictions  on  C,  given  bv  free  center  of  mass  motion  are  then 
k  v 


4T  ckf 

f  u  z 


L  •'m"  (Gkf  Xf  -  Gkf'  xh  =  0  . 

f  9  ox  v  uy  x 


0  =  0 


0  =  0 


(3.6a) 


(3.6b) 


=  +  if  the  eigenvector  of  the  transverse  acoustic  mode  is  q^(+  y)  =  <2rrT.  ( l ,  +  i ,  0)  , 


ind  gives  the  restrictions  on  G: 


('  *  '  ►'in  +•  X  G '  /m  =0, 
nx  j  av 


Cki  ,m  +  :•  Gkv  £  -  0. 

■  IX  .1  -  „  ay  9 


"  =  +  i|i 


(3.8a) 


(3.8b) 


Thus  we  see  that  the  C  symmetry  distinguishes  the  normal  modes  at  0  =  0 

as  symmetric  and  antisymmetric,  at  !l  ^  0  the  C ^  operation  mixes  the  normal  mode 

.it  .ind  to  again  form  svmmetrii-  and  antisymmetric  rush’s.  (See  eqn.  (3.4)). 

It  is  known  that  time  reversal  combined  with  an  inversion  element  can  often  he 

7 

used  to  make  the  dynamic  matrix  and  its  eigenvector  real. 


IV.  Selection  Rules 


The  infrared  and  Raman  selection  rules  for  helical  molecules  have 

hei'ii  discussed  by  lliggs^  and  the  results  are  the  following: 

infrared  absorption  can  only  occur  at  0  =  0,  +  p,  and  Raman  scattering, 

it  i)  =  0,  +  ip,  +  2 ip.  In  general  these  selection  rules  are  applicable 

to  double  helical  molecules.  For  the  molecule  possessing  symmetry 

these  selection  rules  are  affected  by  the  symmetry  and  the  results 

are  Listed  in  Tables  II  and  III  where  x  =  x  +  iy,  x  =  z. 

+  o 

At  a  =  0  the  symmetric  modes  produce  no  y  and  z  components  of 
induced  dipole  and  the  antisymmetric  modes  produce  no  x  component  of 
induced  dipole.  At  b  /  0,  the  normal  mode  at  0  and  -h  are  mixed  by  the 
C  operation.  The  time  reversal  represents t ions  have  been 

classified  by  Lax  as  of  3  types:  real,  pseudoreal  and  complex.  In 
the  vibrational  mode  problem  a  real  representat ion  means  no  double 
degeneracy  due  to  time  reversal,  the  subspace  spanned  by  q(P)  and  the 
subspace  spanned  by  q(-b)  are  identical.  Thus  transitions  that  occur 
at  11  can  also  occur  at  These  results  are  reflected  in  Table  II  and  I1T. 

Furthermore  the  C.j  axis  provides  a  special  direction.  For  Z-DNA  at 
I  lie  pitch  angle  p  r  -t>0"  there  are  b  such  special  directions  in  space. 

When  incident  tight  is  polarized  along,  these  directions  additional 
transition  will  occur  at  11  0.  These  transitions  are  generally  not  found 

lor  a  helical  molecule  without  C.}  symmetry.  The  last  line  of  Tables  II  and  II 
I  i  s  t  s  l  best*  I  rails  i  I  i  mis  . 


T3T3*  W  'JV'.V  A  W  'AW.',1  «.«  v.’  V  TT 
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V.  Conclusions 

Al  lcrfj.it  inn  polymers  of  DNA  have*  .,,,  addirional  symmetry,  a 

i  symmetry  about  the  x-axis.  This  is  in  addition  to  the  screw  axis  or 
helical  symmetry  found  in  all  repeating  DNA’s.  This  additional  symmetry  can 
be  used  along  with  the  helical  symmetry  and  the  ubiquitus  time  reversal  in¬ 
variance  to  greatly  simplify  helix  lattice  calculations  for  these 
alternating  polymers.  In  this  paper  we  have  detailed  the  way  to  make  use  of 
Lhis  simplification.  We  show  how  to  reduce  the  dimension  of  the  dynamical 
matrix  where  possible,  to  greatly  reduce  the  diagonalization  time  needed  for 

solution.  In  practice  we  have  reduced  the  calculation  time  for  the  modes  of 

I  4 

an  alternating  methylated  Z-DNA  by  a  factor  of  ten. 

We  have  also  examined  the  implications  of  the  additional  symmetry  in 
changing  the  infrared  and  Raman  selection  rules.  We  find  that  a  new  absorp¬ 
tion  mode  is  possible  in  these  polymers  with  symmetry.  The  additional 
absorption  mode  should  cause  greater  total  absorption  for  the  alternat¬ 
ing  polymers  than  for  simple  homopolymers.  The  increase  is  strongest  for 
light  polarized  transverse  to  the  helix  axis.  In  crystalline  DNA  where  the 
t>  x-axis  directions  can  be  specified  the  selection  rules  may  he  ver i f ied 
direct lv  by  increased  absorption  at  these  special  directions  traverse  to  the 
helix  axis. 

Acknowl edgement  s 

This  work  has  been  supported  in  part  by  N1H  grant  GM  24443  and  ONR  grant 
NtKHll  4-8D-K-02 52 . 


-12- 

References 

1.  K.  C.  I.ii,  E.  W.  Prohofsky  and  L.  I..  Van  Zandt,  Biopolymers  16,  2491  — 

2r>06  (1977). 

2.  K.  V.  Devi-Prasad  and  E.  W.  Prohofsky,  J.  Biomolec.  Struct,  and  Dynamics 
2,  627-635,  (1984). 

3.  J.  C.  Decins  and  R.  M.  Hexter,  "Molecular  Vibrations  in  Crystal", 
(McGraw-Hill,  New  York,  1977). 

4.  C.  Y.  I.iang,  J.  Mole,  Spectr.  1^,  61  (1957). 

5.  E.  W.  Prohofsky,  "Motional  Dynamics  of  the  DNA  Double  Helix,"  Biomolecular 
Stereodynamics  IV,  ed.  by  R.  H.  Sarmar  and  M.  H.  Sarmar,  (Adenine  Press, 
New  York,  1986),  and  references  therein. 

6.  J.  M.  Evster  and  E.  W.  Prohofsky,  Biopolymers  2505  (1974). 

7.  M.  Lax,  "Symmetry  Principles  in  Solid  State  and  Molecular  Physics", 

(Wiley,  New  York,  1974). 

8.  .1.  E.  Rosenthal  and  G.  M.  Murphy,  Rev.  Mod.  Phys.  8,  317  (1936). 

9.  1).  A.  Pearlman  and  S.  H.  Kim,  Biopolymers  U_,  327  (1985). 

10.  W.  Saenger,  "Principles  of  Nucleic  Acid  Structure",  (Springer-Verlag, 

New  York,  1984). 

11.  Arnott  and  Chandrasekaran ,  private  communication  (1985). 

12.  Arnott  and  Chandrasekaran,  private  communication  (1987). 

13.  P.  W.  Higgs,  Proc.  Rov.  Soc.  A 2 2 0 ,  472  (1953). 

14.  X.  fl.  Hun  and  E.  W.  Prohofsky,  "Normal  Mode  Calculation  for  Methylated 

5  5 

/-DNA  pnlv(dC-m  dC)  (dG-m  dC)",  to  be  published  in  Biopolymers. 


x,y  real  and 
arbitrary 


